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Abstract 

A graph (digraph) G = iy, E) with a set T C 1/ of terminals is caUed 
inner Eulerian if each nonterminal node v has even degree (resp. the numbers 
of edges entering and leaving v are equal). Cherkassky JT] and Lovasz 
showed that the maximum number of pairwise edge-disjoint T-paths in an inner 
Eulerian graph G is equal to ^ X^seT where A(,s) is the minimum number 
of edges whose removal disconnects s and T — {s}. A similar relation for inner 
Eulerian digraphs was established by Lomonosov Considering undirected 
and directed networks with "inner Eulerian" edge capacities, Ibaraki, Karzanov, 
and Nagamochi JOl showed that the problem of finding a maximum integer 
multiflow (where partial flows connect arbitrary pairs of distinct terminals) 
is reduced to O(logT) maximum flow computations and to a number of flow 
decomp osit ions . 

In this paper we extend the above max-min relation to inner Eulerian bidi- 
rected and skew- symmetric graphs and develop an algorithm of complexity 
0(yii^ log Tlog(2 + V'^/E)) for the corresponding capacitated cases. In par- 
ticular, this improves the bound in JHl for digraphs. Our algorithm uses a fast 
procedure for decomposing a flow with 0(1) sources and sinks in a digraph into 
the sum of one-source-one-sink flows. 
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1 Introduction 



A graph (digraph) G = (V, E) with a distinguished subset T of nodes is said to be 
inner Eulerian if each node v E V — T has even degree (resp. the indegree and 
outdegree of v are equal). The nodes in T and m V — T are called terminals and 
mner nodes, respectively. A simple path in G is called a T-path if its ends are distinct 
terminals and the other nodes are inner. There is a nice max-min relation established 
by Cherkassky pP and Lovasz ^2] for graphs, and by Lomonosov ^3] for digraphs: 

(1) if (G, T) is inner Eulerian, then the maximum number ug,t of pairwise edge- 
disjoint T-paths is equal to | YIs^t ^g,t{.s). 

Here and later on, for a subset X of nodes, 5{X) = 6g{X) denotes the set of edges 
with one end in X and the other in ^ — X, called the cut induced by X. For s E T, 
we refer to a subset X G V with X (IT = {s} as an s-set. Then Xg,t{s) is defined to 
be the minimum cardinality |5(X)| among the s-sets X. 

The above max-min relation has an obvious extension to the capacitated case. 
Given a nonnegative integer function c : E ^ Z+ of edge capacities, let us say that 
the triple (network) [G, T, c) is inner Eulerian if for each inner node v, the total 
capacity of edges incident with v is even when G is a graph, and the total capacity 
of edges entering v is equal to that of edges leaving v when G is a digraph. Then ((H) 
yields the following relation for an inner Eulerian [G, T, c): 



where Ac,t(s) denotes the minimum cut capacity c{6{X)) among the s-sets X, the 
maximum is taken over all collections of T-paths Pi, . . . ,Pk along with nonnegative 
integer weights ai, . . . ,ak G 1^+ that satisfy the packing condition 



and val(jF) denotes the total value of JF, defined to be ai + . . . + a^^. (Hereinafter 
for a function g : S ^ M. and a subset S' C 5*, g{S') stands for J2eeS' di^)-) 

A collection JF consisting of T-paths Pi with real weights G M+ that obeys 
is called a free multiflow (the adjective "free" is used to emphasize that any pair of 
distinct terminals is allowed to be connected by a path, i.e., the commodity graph in 
the multiflow maximization problem is complete). A multiflow achieving the equality 
in Q is called maximum. Thus, whenever {G, T, c) is inner Eulerian, there exists an 
integer maximum free multiflow (i.e., having the weights of all paths integral). 

Cherkassky (TJ showed that such a multiflow in an inner Eulerian undirected net- 
work can be found in strongly polynomial time. Subsequently much faster algo- 
rithms both for graphs and digraphs have been developed. They apply a "divide-and- 
conquer" approach in which a current network {G,T,c) with |T| > 4 is recursively 

^Originally relation ^ was stated for fractional multiflows T in ^2]) with a flaw in the proof. 
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for all e G E, 
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replaced by two networks (G", T', c') and (G", T", c") such that |T'|, |T"| < [|T|/2] +1. 
Originally, such an approach was applied in [TT] to find, in 0{(j){V,E) logT) time, a 
half-integer maximum multifiow in a graph G with integer edge capacities (but not 
guaranteeing integrality in the inner Eulerian case). Hereinafter, in notation involving 
functions of numerical arguments or time bounds, we indicate sets for their cardinal- 
ities, and 0(n, m) stands for the complexity of an algorithm for finding a maximum 
fiow in a network with n nodes and m edges. 

This algorithm was improved and extended in PU] so as to find an integer max- 
imum free multifiow in an inner Eulerian undirected network in the same time 
0{(j){V, E) logT), and in an inner Eulerian directed network in 0{(j){V, E) logT+V'^E) 
time. 

Remark 1. The inner Eulerianness condition is important. Withdrawing it makes 
the undirected problem more difficult, though still polynomially solvable in the non- 
capacitated case (a max-min relation is due to Mader and an original polynomial 
algorithm is due to Lovasz pUSj), and makes the directed noncapacitated problem 
NP-hard already for two terminals |3]. 

The purpose of this paper is to extend the above theoretical and algorithmic results 
to bidirected graphs. (This sort of nonstandard graphs was introduced by Edmonds 
and Johnson j2] in connection with one important class of integer linear programs 
generalizing problems on fiows and matchings; for a survey, see also [1^1 llHj.) 

Recall that in a bidirected graph G = {V, E) three types of edges are allowed: (i) 
a usual directed edge, or an arc, that leaves one node and enters another one; (ii) 
an edge from both of its ends; or (iii) an edge to both of its ends. When u = v, the 
edge becomes a loop; in what follows we admit only loops of types (ii) and (iii) (as 
loops of type (i) do not affect our problem and can be excluded from consideration). 
A nonloop edge entering a node v contributes one unit to the indegree deg™(t') of v, 
while a loop of type (iii) at v contributes two units to deg"(f ); the outdegree deg°"*(f ) 
of V is specified in a similar way. Edges e, e' connecting nodes u, v are called parallel 
if e enters u if and only if e' does so, and similarly for v. If G has no parallel edges, 
then \E\ < 2\V\'^. An instance of bidirected graphs is drawn in Fig. ^ 

The notion of inner Eulerianness for a bidirected graph G with a set T of terminals 
is analogous to that for usual digraphs: deg™(f) = deg°"*(f) for all inner nodes v. 
Inner Eulerian triples {G,T,c), where c : E Z+, are those that turn into inner 
Eulerian pairs (G', T) when each edge e is replaced by c(e) parallel edges. 

In order to be able to extend the above results to bidirected graphs, we need to 
admit T-paths with restricted self-intersections. (In undirected or directed graphs, 
when one node is reachable from another one by a path, then it is reachable by a simple 
path as well, but this need not hold in a bidirected graph.) A walk in a bidirected 
graph G is an alternating sequence P = (s = fo, Ci, f i, . . . , e^, f ^ = t) of nodes and 
edges such that each edge Cj connects nodes and Vi, and for z = 1, . . . , — 1, the 
edges Ci, Cj+i form a transit pair at Vi, which means that one of Cj, Cj+i enters and 
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Figure 1: A bidirected graph instance 



the other leaves Vi. Note that ei may enter s and Cfc may leave t; nevertheless, we 
refer to P as a walk from s to t, or an s-t walk. P is a cycle if Vq = Vk and the pair 
ei, Cfc is transit at fo; a cycle is usually considered up to cyclically shifting. Observe 
that an s-s walk is not necessarily a cycle. By a path in a bidirected graph G we 
mean an edge-simple walk P, i.e., a walk with all edges different. Similar to usual 
graphs/digraphs, a T-path (a T-walk) is meant to be a path (resp. walk) whose ends 
are distinct terminals and the other nodes are inner. 

Define the number \g,t{s) for s G T as before. We show the following 

Theorem 1.1 Property (OJj remains valid for a bidirected graph G and the set of 
T-paths in it. 

Remark 2. In this theorem it suffices to consider only minimal T-paths, where a 
path (edge-simple cycle) P = {vq, Ci, Vi, . . . , e^, Vk) is called minimal if no part of P 
from Vi to Vj with i < j (resp. < j — i < k) forms a cycle. (A minimal path/cycle 
need not be simple but it passes any node of G at most twice.) Moreover, one can 
consider only those T-paths whose induced bidirected graphs contain no cycle. (In 
the underlying undirected graph of such a path, each edge belongs to at most one 
circuit.) 

A usual digraph is a special case of bidirected graphs and Theorem II. II generalizes 
the above-mentioned result in jT^. Also there is a natural correspondence between 
the T-paths in an undirected graph G and the minimal T-paths in the bidirected 
graph G' formed from G as follows: direct each edge of G from both of its ends, and 
for each inner node v, assign [deg(f)/2] loops entering (twice) v. Then a T-path 
P in G is turned into a T-path in G' by adding one loop to each intermediate node 
of P. Moreover, (G", T) is inner Eulerian if (G, T) is such. Due to this correspon- 
dence. Theorem 11.11 generalizes the above-mentioned Cherkassky-Lovasz' result for 
undirected graphs as well. 

Like the pure graph and digraph cases, one can reformulate Theorem II. II in capac- 
itated terms: relation ^ concerning integer free multifiows remains valid when G 
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is bidirected and {G, T, c) is inner Eulerian. In this case one should consider T-walks, 
rather than T-paths, and refine the packing condition Q as 

k 

(4) ^ai^i(e) < c(e) for all e G -E, 

i=l 

where ni{e) is the number of occurrences of an edge e in a walk Pj. Thus, the above 
problem for undirected and directed networks is generalized as: 

(P) Given an inner Eulerian network {G,T,c), where G is bidirected, find a collec- 
tion (free multiflow) T of T-walks Pi, . . . ,Pk with weights ai, . . . ,ak G that 
satisfies ^ and maximizes the value val(jF) := ai + . . . + a^- 

Remark 3. Let X be an arbitrary subset of nodes of G. One can modify G as follows: 
for each node f G X and each edge e incident with v, reverse the direction of e at v. 
Also for an arbitrary arc e incident with a terminal s, one can reverse the direction of 
e at s. Both transformations preserve the inner Euleriannes of {G, T, c) and the set 
of T-walks. Therefore, problem (P) remains equivalent under such transformations. 

Another appealing class of nonstandard graphs was introduced by Tutte [20] who 
originated a mini-theory, parallel to 2j in a sense, aiming to unify and generalize 
flow and matching problems. These are so-called skew- symmetric graphs (or anti- 
symmetrical digraphs, in Tutte's terminology), digraphs with involutions on the nodes 
and on the arcs which reverse the orientation of each arc (a precise definition is given in 
Sectional). His and other researchers' study of structural and optimization problems 
on skew-symmetric graphs has resulted in a number of interesting theorems, methods 
and applications. 

There is a close relationship between skew-symmetric and bidirected graphs, and 
typically results on the former can be reformulated for the latter, and vice versa. So is 
for the problem of our study, too. We take advantage from both representations. The 
language of bidirected graphs is more preferable for us to work in the non-capacitated 
case; we prove Theorem 11.11 direct Iv and obtain its analog for skew-symmetric graphs 
as a corollary. On the other hand, we prefer to deal with skew-symmetric graphs 
in algorithmic design for the capacitated case. (A serious reason is that a fiow in a 
bidirected network is defined as a packing of T-paths and we do not see reasonable 
alternative settings for it, while a fiow in a skew-symmetric network can be given in a 
more compact form, via a function on the arc set.) Some facts about skew-symmetric 
fiows and technical tools elaborated for them help us to devise a fast algorithm for 
the skew-symmetric analog of problem (P) concerning integer skew- symmetric free 
multiflows in an inner Eulerian skew-symmetric network. This yields a fast algorithm 
for (P) as well. 

The paper is organized as follows. In Section |21 we prove Theorem 11.11 (which 
is relatively simple) relying on the fact that an inner Eulerian bidirected graph can 
be decomposed into cycles and paths with both ends in T. Section [HI explains the 



4 



correspondence between bidirected and skew-symmetric graphs, reviews some known 
results about the latter (in particular, Tutte's result on symmetric decompositions of 
skew-symmetric flows) and gives a skew-symmetric analog of Theorem ll.il Section E] 
develops an algorithm for flnding a maximum integer skew- symmetric free multiflow in 
an inner Eulerian skew-symmetric network. It relies on a general approach in JT] and 
some ingredients from and attracts additional combinatorial ideas and techniques. 
As a consequence, problem (P) is solved in time 0(yi^^ log T log(2 + V'^/E)) (if the 
0(nm log(2 + n^/m) )-algorithm of Goldberg and Tarjan [0] is applied for flnding a 
maximum flow in a directed network with n nodes and m arcs). This improves the 
bound in for digraphs. To achieve this bound, we use a faster procedure for the 
particular flow decomposition problem: given an integer flow / with 0(1) sources 
and sinks in a digraph with n nodes and m arcs, decompose / into the sum of integer 
flows, each connecting one source to one sink. The procedure developed in Sectional 
solves this problem in 0(m log(2 + n^/m)) time. In the concluding Section this 
procedure is extended to symmetric flows in skew-symmetric graphs (it is not used in 
the algorithm for (P) but may be of interest for other applications). 

2 Proof of Theorem 11.11 

Let G = (y, E) be an inner Eulerian bidirected graph with a terminal set T. One may 
assume that G has no loops incident with terminals. Since (G, T) is inner Eulerian, 
for each inner node one can choose a set 7r(f) of transit pairs at v so that each 
non-loop edge incident with v occurs in exactly one pair and each loop at v (if any) 
occurs in two pairs. The collection {t^{v) : v &V — T} determines a decomposition 
of (the edge-set of) G into a collection C of edge-simple cycles and a collection V 
of paths with both ends at T. More precisely, each edge e E E belongs to exactly 
one member P of C U V and satisfles the following condition: for each end v of e, 
if f G V — T and {e, e'} G vr(f), then either e,v,e' or e',v,e are three consecutive 
elements in P, while if u G T, then P begins with v,e or ends with e,v. Note that 
all nodes of any cycle in C and all intermediate nodes of any path in V are inner. So 
each path in P is a T-path unless it connects equal terminals. When needed, we may 
reverse some paths in V. 

For s G T, let Vg (Qs) denote the set of paths in V with exactly one end (resp. 
with both ends) at s. Since \Vs\ + 2|Qs| = deg(s) > Xg,t{s) (where deg(f) is the full 
degree deg°"*(t') -|- deg™(f) of v), the theorem becomes trivial when all sets Qs are 
empty. In a general case, we try to transform the decomposition so as to increase the 
"useful value" ri(V,C) := J2s£T by applying a certain augmenting approach. 

Consider s G T and assume, w.l.o.g., that all paths in Vs begin at s. Let L = 
(xo, Xi, . . . , Xq) be a sequence of distinct nodes such that 

(5) either L = {s}, or xq belongs to a path in Qg, and for i = 1, . . . , g, the nodes 
Xi, Xi-i occur in a cycle in C or occur in this order in a path in Vg- 
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We say that L is augmenting if Xq belongs to a path in V having both ends in T — {s}. 
Consider two cases. 

Case 1. There is no augmenting sequence for s. Let Xg be the set of all nodes 
occurring in sequences L as in (0). Clearly fl T = {s}. Consider an edge e of the 
cut 5{Xs)] let Mjf be the ends of e in Xs and V — X^, respectively. Observe that e 
belongs to neither a cycle in C nor a path ixiV — Vs- Also (by ©) if e belongs to a 
path P & Vs, then all nodes of P from s to (the last occurrence of) u are contained 
in Xg, i.e., P traverses the cut S{Xs) exactly once. This implies \Vs\ = \S{Xs)\. 

Hence, if none of terminals admits an augmenting sequence as above, then the 
number of T-paths in V is at least ^ J2seT ^g,t{s), as required. 

Case 2. An augmenting sequence L = {xq, xi, . . . , Xq) for s exists. Let L be chosen 
so that no proper subsequence in it is augmenting. Then: 

(6) any cycle in C meets at most two nodes in L and these nodes are consecutive 
in L; 

(7) if a path P E Vs contains a node Xi, then the part of P from s to (the last 
occurrence of) Xi can contain at most one node Xj with j > i; moreover, if such 
an Xj exists then j = i + 1. 

We transform (V, C) along L, step by step, as follows. Choose a. Q E Qs containing 
Xq. At the first step, if (a) xo,xi belong to a cycle C G C, then we combine Q and 
C into one s-s path. And if (b) Xo,Xi belong to a path P E Vs from s to t, say, 
and if Xi occurs in P earlier than xq, then we replace Q by the concatenation of the 
part P' of P from s to (the last occurrence of) xq and the part Q' of Q from xq to 
s, and replace P by the concatenation of the rest of Q (from s to xo) and the rest 
of P (from Xq to t). (We assume, w.l.o.g., that the last edge of P' and the first edge 
of Q' form a transit pair at Xq; otherwise reverse Q.) As a result, we obtain an s-s 
path, denoted by Q as before, that contains Xi. In case (a), the cycle C vanishes, 
and in case (b), the new path P goes from s to t as before, and its part from xq to 
t preserves. This together with (jHl) and ((Zj) implies validity of ^ for the remaining 
sequence (xi, . . . ,Xq); moreover, (jH)) and ((Zj) are maintained as well. At the second 
step, we consider the pair Xi,X2 and act in a similar way, and so on. 

Eventually, after q steps, the current s-s path Q contains the node Xg. Since L is 
augmenting, Xq also belongs to some t-p path R E V with t,p E T — {s} (possibly 
p = t). Now splitting Q, R at Xq and concatenating the arising four pieces in another 
way, we obtain two T-paths, one connecting s and p and the other connecting s and 
t. This gives a new decomposition {V, C) of G having a larger value of r], and the 
theorem follows. ■ 

In fact, the above proof is constructive and prompts a polynomial algorithm for 
finding a maximum number of pairwise edge-disjoint T-paths in an inner Eulerian 
bidirected graph. A more efficient and more general algorithm (dealing with the 
capacitated case) is described in Section HI 
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3 Skew-Symmetric Graphs 



This section contains terminology and some basic facts concerning skew-symmetric 
graphs and explains the correspondence between these and bidirected graphs. For a 
more detailed survey on skew-symmetric graphs, see, e.g., 171 15]. 

A skew- symmetric graph is a digraph G = {V, E) endowed with two bijections 
ay, (Je such that: ay is an involution on the nodes (i.e., av{v) ^ v and av^cviv)) = v 
for each v ^ V), as is an involution on the arcs, and for each arc e from u to v, o"£;(e) 
is an arc from o"y(f ) to av{u). For brevity, we combine the mappings ay, cte into one 
mapping a on V U E and call a the symmetry (rather than skew-symmetry) of G. 
For a node (arc) x, its symmetric node (arc) a{x) is also called the mate of x, and 
we will often use notation with primes for mates, denoting cr{x) by x'. Obviously, 
deg'^'iv) = deg°''\v') for each veV. 

We admit parallel arcs, but not loops, in G. Observe that if G contains an arc e 
from a node v to its mate v', then e' is also an arc from v to v' (so the number of arcs 
of G from v to v' is even and these parallel arcs are partitioned into pairs of mates). 

By a path (circuit) in G we mean a simple directed path (cycle), unless explicitly 
stated otherwise. The symmetry a is extended in a natural way to paths, subgraphs, 
and other objects in G. In particular, two paths or circuits are symmetric to each 
other if the elements of one of them are symmetric to those of the other and go in 
the reverse order: for a path (circuit) P = {vq, Ci, Vi, . . . , e^, Vk), the symmetric path 
(circuit) cr(P) is {v[, e'l^, v'i._i, . . . , e[, v'q) . One easily shows that G cannot contain 
self- symmetric circuits (cf. jH]). 

Following terminology in , a path or circuit in G is called regular if it contains 
no pair of symmetric arcs (while symmetric nodes in it are allowed) . For a function 
h on E, its symmetric function h' is defined by h'{e') := h{e), e E E, and h is called 
(self-) symmetric if h = h' . 

For a function f : E ^ and a node v E V, define 

divfiv) := 5^(/(e) : e G 6°-\v)) - ^^(/(e) : e G 6'-{v)), 

(the divergency of / at f ), where 6°^^{v) (5™(f)) denotes the set of arcs of G leaving 
(resp. entering) v. Let / be nonnegative, integer-valued and symmetric, and let S be 
a subset of nodes not intersecting S' = <7{S). When divj(t>) is nonnegative at each 
V E S and zero at each v E V — {S U S'), f is said to be an IS- flow (integer symmetric 
flow) from S to S' . The value val(/) of / is ^sgsdiv/(s). By a multiterminal version 
of a theorem due to Tutte |20j, an IS-fiow / from S to S' has an integer symmetric 
decomposition. This means that 

(8) / is representable as / = aix^^ + ctix^^ + . . . + akX^*' + oikX^'^i where for 
i = 1, . . . , A;, Pj is a path from S to S' or a circuit, P/ is the path (also going 
from S to S') or circuit symmetric to Pj, and G Z+. 
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Here denotes the incidence vector of the arc-set of a path/circuit P, i.e., for e G -E, 
X^(e) = 1 if e belongs to P, and otherwise. Note that paths/circuits in (jH)) need 
not be regular. Considering as the weight of Pi and of P/, observe that the total 
weight of paths from S to S' is equal to val(/). Similar to flow decomposition in 
usual digraphs, an integer symmetric decomposition of an IS-fiow / can be found in 
0{VE) time. 

Let {G, T) be inner Eulerian, where the terminal set T is (self-)symmetric. Take 
a partition {S, S' = <j{S)} of T such that deg°"*(s) > deg™(s) for all s G S. Since the 
all-unit function f on E represents an IS-slow from S to S', (jH)) implies that 

(9) there exists a symmetric collection V of circuits and paths from S to S' in 
G such that the members of V are pairwise arc-disjoint and cover E, and each 
terminal s G S" is the beginning of exactly deg°"*(s) - deg'''(s) paths in V. 

Moreover, the members of V are regular (for if some P G P contains mates e, e' G E, 
then e, e' are in cr(P) as well, which is impossible). 

Next we explain the correspondence between skew-symmetric and bidirected 
graphs (cf. 0| Sec. 2]). For sets X,A,B, we may use notation X = A U B when 
X = AU B and A fl P = 0. Given a skew-symmetric graph G = (V, E), choose an 
arbitrary partition vr = {Vi,V2} of V such that V2 is symmetric to Vi. Then G, vr 
determine bidirected graph H with node set Vi whose edges correspond to the pairs 
of symmetric arcs in G. More precisely, arc mates a, a' of G generate one edge e oi H 
connecting nodes u,v E Vi such that: (i) e goes from m to f if one of a, a' goes from u 
to V (and the other goes from v' to u' in V2); (ii) e leaves both u, v if one of a, a' goes 
from u to v' (and the other from v to u'); (iii) e enters both u, v if one of a, a' goes 
from m' to V (and the other from v' to u) . In particular, e is a loop if a, a' connect a 
pair of symmetric nodes. 

Conversely, a bidirected graph H with node set Vi, say, determines skew- 
symmetric graph G = (V, E) with symmetry a as follows. Take a copy (T(f ) of 
each element v of Vi, forming the sets V2 := {<j{v) : f G Vi} and V := ViU V2. For 
each edge e oi H connecting nodes u and v, assign two "symmetric" arcs a, a' in G so 
as to satisfy (i)-(iii) above (where u' = a{u) and v' = cr{v)). An example is depicted 
in Fig. H 

Remark 4. A bidirected graph generates one skew-symmetric graph, while a skew- 
symmetric graph generates a number of bidirected ones, depending on the partition 
TT of that we choose in the first construction. The latter bidirected graphs are 
produced from each other by the edge reversing transformation with respect to a 
subset of nodes as indicated in Remark 3 in the Introduction, so they are equivalent 
for us. 

A terminal set S in H generates the symmetric terminal set T := S U cr(5') in G, 
and vice versa. One easily checks that {H, S) is inner Eulerian if and only if (G, T) is 
such. Also there is a correspondence between the S'-paths in H and certain T-paths 
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Figure 2: Related bidirected and skew-symmetric graphs 



in G. More precisely, let r be the natural mapping of U to Vi U E{H) (where 
E{H) is the edge set of H). Each walk (cycle) P = {vq, ai,vi, . . . , a^, Vk) in G induces 
the sequence r(P) := {t{vo), T(ai), t{vi), . . . , r(afc), T{vk)) of nodes and edges in H. 

Conversely, for a walk (cycle) Q = {wq, ei,Wi, . . . , e^, Wk) in H, form the sequence 
t{Q) := {vq, ai, f 1, . . . , Ofc, Vk) of nodes and arcs in G by the following rule: 

(R) Vq := Wq if Ci Icavcs Wq, and Vq := o-{wo) if Ci enters Wq; and for i = 1, . . . , fc: 
(a) if Cj leaves then Oj is the arc in T~^{ei) that leaves and fj is the 
head of af, (b) if Cj enters then Oj is the arc in r~^(ej) that leaves cT(wj_i), 
and Vi is the head of Oj. 

(When is a loop, the arcs in r~^(ei) are parallel, and the arc in this set is chosen 
arbitrarily.) It is not difficult to conclude that (R) provides: 

(10) for a walk (cycle) Q in H, 

(i) t{Q) is a walk (cycle) in G and t{t{Q)) = Q; 

(ii) if Q is edge-simple and minimal (see Remark 2 in the Introduction), then 
t{Q) is a regular path (circuit). 

Also the walk (cycle) reverse to Q determines the walk (cycle) in G symmetric to 
t{Q) (up to the choice of arcs Oj for loops Cj). The corresponding converse properties 
to those in (fTUj) also take place. 

Let us say that a T-walk P from s to t in G is essential if t is different from 
ct(s). Thus, we have a natural bijection between the essential regular T-paths in G 
(considered up to parallel arc mates) and the minimal S'-paths in H. This gives 

(11) I'CT = '2l^H,S, 

where T^g,t is the maximum cardinality of a symmetric collection of pairwise arc- 
disjoint essential T-paths in G. 
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Note also that for a terminal s E S and an s-set X in H, each edge of the cut 
6h{X) generates two arc mates in the symmetric cut Sg{X L\a{X)) in G. Therefore, 

(12) Ag,t(s) = 2\h,s{s) for each s E S, 

where Xg.t{s) is the minimum cardinality of a symmetric cut in G separating {s, s'} 
and T — {s, s'}. 

In view of relations (jlip and (jl2p . Theorem II. II is equivalent to the following 

Corollary 3.1 For a skew-symmetric graph G = {V,E) with a symmetric set T of 
terminals, if {G,T) is inner Eulerian, then ug,t = jYlses ^g,t{s) . 

In the capacitated case, we are given a symmetric function c : E Z_|_ of arc 
capacities in a skew-symmetric graph G = (V, E) with a symmetric set T = 5 U 5" 
of terminals. By an integer symmetric free multiflow (or, briefly, an IS-multifiow) in 
the network {G, T, c) we mean a collection F of integer flows fst for the ordered pairs 
(s, t) of distinct terminals in S such that: (a) / = fst is a flow from {s, s'} to {t,t'} 
(i.e., divf{v) is nonnegative for v = s,s', nonpositive for v = t,t', and otherwise); 
(b) each fst is symmetric to fts] and (c) F is c- admissible, i.e., 

^ /st(e) < c(e) for each e E E. 

The (total) value val(F) of F is Ylst^^Kfst)- The problem is: 

(PS) Given an inner Eulerian network {G,T,c), where G is a skew- symmetric graph 
and T and c are symmetric, find a maximum IS-multiflow, i.e., an IS-multiflow 
F maximizing val{F) . 

To see how this problem is related to (P), let (G, T) correspond to {H, S), where H 
is bidirected. Let c be the corresponding capacity function in H, i.e., c(e) = c(a) for 
an edge e E E{H) and its images a, a' in G. The inner Eulerianness of {G, T, c) implies 
that of {H, S", c), and vice versa. Given an IS-multifiow F in {G, T, c), represent each 
flow fst in the path packing form: 

(13) fst = ctix^^ + . . . + a/cX^*! where = a{Pi) E and Pi is a circuit or a 
(simple) path from {s, s'} to {t,t'}. 

We assume that the representation of each flow fst is symmetric to that of fts- Then 
the set V of (essential) T-paths in these representations is symmetric, with a{P) = 
<y{P') for each P E V, and we have val(F) = "^{(yiP) '■ P E V). Now each pair 
P,P'eV of path mates determines an 5- walk P in H (considered up to reversing), 
and taking together these paths P with weights a(-P), we obtain a multiflow JF in 
{H,S,c) satisfying val(J-') = |val(F). 

Conversely, let JF = (V, a) be an integer multiflow in {H, S,c), where V consists of 
S'-walks. One may assume that for each edge e of H, no path P eV traverses e twice 
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in the same direction (for otherwise one can remove a cycle from P). Then each s-t 
walk P determines an arc-simple directed walk P from {s,s'} to {t,t'} and its mate 
P' from {t,t'} to {s, s'} in G. Assign a{P) := a{P') := a{P). Let fst be the sum of 
functions ax^ over the obtained walks P from {s, s'} to {t, t'}. Then fts is symmetric 
to /sj. These flows form an IS-multiflow F in {G,T,c) satisfying val(F) = 2val(jF). 

Thus, problems (PS) and (P) (regarding H, S,c) are reduced to each other. In the 
next section we devise an efficient algorithm for finding an optimal solution to (PS) 
and then explain that it can be transformed into an optimal solution to the corre- 
sponding instance of (P) without increasing the time bound. 

4 Algorithm 

In this section we describe an algorithm to solve problem (PS) and estimate its 
complexity. We use terminology and facts from the previous section. 

Let {G = {V, E),T,c) be an inner Eulerian skew-symmetric network. As before, 
we represent the terminal set T as SUS' and associate with [G, T, c) the corresponding 
bidirected network {H,S,'c). One may assume that no arc in G connects a pair of 
terminal mates. Also if G has an arc entering a terminal s G 5*, then replacing its head 
s by s' and symmetrically replacing the tail s' of the symmetric arc e' by s does not 
affect the problem in essence. So we may assume that deg™(s) = for each terminal 
s & S in G. Then any fiow from {s, s'} to {t,t'}, where s,t E S, is essentially a fiow 
from s to t', and its symmetric fiow is a fiow from t to s'; this property will simplify 
technical details in our construction. In terms of H, the latter assumption says that 
each edge incident with a terminal s in if leaves s (cf. Remark 3 in Section ^). 

The algorithm uses a recursion analogous to that in [TT], and the case |S'| = 3 is 
the base in it. We first consider this special case (which generalizes the case \S\ = 2). 

4.1 Case l^l = 3. 

The algorithm for this case uses one auxiliary skew-symmetric graph Gi = {V, Ei). It 
is obtained from G by adding, for each pair v,v' of inner node mates, four auxiliary 
arcs connecting v and v': two arc mates a^,, a'^ going from v to v' and two arc mates 
a„/, a'^, from v' to f , regardless of the existence of such arcs in G. This Gi corresponds 
to the bidirected graph Hi obtained from H by adding two auxiliary loops at each 
inner node f , one leaving v (twice) and the other entering v. 
The algorithm consists of three stages. Let S = {si, S2, S3}. 

At Stage 1, we apply the algorithm for inner Eulerian graphs from ^lOj to find 
a maximum integer free multifiow in the underlying undirected graph H for H hav- 
ing the same set S of terminals and the same capacities a It runs in 0{(f){V, E)) 
time (since \S\ = 0(1)) and outputs (simple) S'-paths Pi, . . . , P^ in H and weights 
tti, . . . , ttfc G Z_|_ satisfying the packing condition w.r.t. a (Recall that is a bound 
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Pi 

Figure 3: A path Pi and its images in the graphs Hi and Gi 

for the apphed max flow algorithm; we assume = Q{VE)). It also outputs pairwise 
disjoint Sj-sets Xj, i = 1,2, 3, such that for each i, the sum of weights aj of paths Pj 
connecting Sj and S — {sj} is equal to c{6H{Xi)). However, some pairs of consecutive 
edges in Pj may be non-transit in H, i.e., Pj is not necessarily a path in H. 

At Stage 2, we transform Pi, . . . , Pk into paths in the auxiliary bidirected graph 
Hi. More precisely, for each Pi = (t>o, ei, vi, . . . , Cq, Vq) and for each non-transit pair 
Cj, Cj+i in it, if both edges Cj, Cj+i enter (leave) Vj, then the element Vj of Pi is replaced 
by the string Vj,i,Vj, where i is the auxiliary loop leaving (resp. entering) Vj. This 
results in minimal S'-paths Pi, . . . ,Pk in Hi. 

Each path Pj and its reverse one are then lifted to Gi (by the method explained in 
Sectional), giving regular T-paths Pi, P/ symmetric to each other. (FigureEliUustrates 
paths Pi, Pi, Pi.) For each pair Sj,Sp {j 7^ p), the functions aiX^' or a^x^^ for the 
paths from Sj to Sp are added up, forming Sj-Sp flow gjp. This gives a symmetric 
collection of six integer flows in Gi, see Figure 0] The c-admissibility of the above 
multiflow in H and the fact that each path Pj is regular imply that the total flow 
though each arc e of G does not exceed c(e). Also the fact that the cuts 5//(Xj) are 
saturated implies that 

(14) for i = 1,2, 3, the arcs in (5°'^*(Xj) are saturated by 5'j,j--i + gi,i+i, and sym- 
metrically, the arcs in (5™(Xj) are saturated by 5'j_i,j + gi+i,i, 

where Xj := Xj U X- and the indices are taken modulo 3. So the IS-multiflow 
consisting of these six flows has maximum value. 

At Stage 3, we improve the above flows gij in Gi by reducing their values on 
the auxiliary arcs to zero, eventually obtaining the desired multiflow in G. In view 
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Figure 4: Flows gij in Gi 

of (dH), for i = 1,2,3, one may assume that (^j.i^j+i and gi+i^i-i take zero values on 
all arcs of the subgraph (Xj) of Gi induced by Xi. 

Take the residual capacities A(e) := c{e) — '^-j gij{e) of arcs e E E. The divergency 
of c (w.r.t. E) and of each Qij (w.r.t. Ei) at any inner node is zero, therefore, 

(15) divA(f ) = / gijiv, v') for each v E V — T, 

where g{v,v') denotes g{ay) + g{a'^) — g{a^') — g{a'^,) (recall that the auxiliary 

arcs from v to v'). The function A on is nonnegative, integer- valued and symmetric. 
Also JIHl) and 

(16) gij{v, v') = —gij{v', v) = gji{v, v') for each v eV — T 

imply that divA(f ) is even for each v E V—T. Hence we can extend A to the auxiliary 
arcs so as to obtain an IS-flow in {Gi,T). (The extended A satisfies A{v,v') + 
dijiv, v') = for each v eV - T.) 

Notice that A(e) = for each arc e in the cut S{Xi), i = 1,2, 3, by (fT^ . Therefore, 
the restriction Aj of A to the set Ai of arcs of the subgraph (Xi) is an IS-flow from 
Si to In its turn, the restriction Aq of A to the set Aq of arcs with both ends in 
W := V^ — (Xi UX2 UX3) is an integer symmetric circulation in the subgraph {W, Aq). 
(Recall that the sets Xi,X2,X3 are pairwise disjoint.) 

We start with getting rid of nonzero arc values of the above flows on the auxiliary 
arcs within the subgraph (Xi). To this aim, apply the integer symmetric decompo- 
sition procedure to Ai (cf. (jHl)) to represent it as the sum of integer si-s'^-flows h, h', 
where h' is symmetric to h. Combine g := gi2 + gis + h and g' := g2i+ g^i + h' (where 
h, h' are formally extended by zeros on Ei — Ai). Then g is an integer flow from Si 
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to S', and g' is the flow from S to s[ symmetric to g. Also 

g{v, v') = gi2{v, v') + gis{v, v') + h{v, v') = for each v E Xi - {si, s[}, 

in view of h{v,v') = |Ai(w,w'), A{v,v') + Y.ij 9ij{v,v') = and 5(23(6) = 532(6) = 
for all e E Ax. So we can reduce 5, g' to zero on all auxiliary arcs in (Xi). Now using 
standard flow decomposition, we represent the new flow g as the sum of three integer 
flows /i,/2,/3, from s\ to s^, from si to s^-, and from s\ to S3, respectively. Note 
that g{e) = for each e G (5™(Xi) implies that /i is zero on all arcs of the cut b{X\). 
Update 5i2 := /2 and 513 := /3; the flows 521 and 531 are updated symmetrically. 
Then the resulting four flows together with the remaining flows 523, ^732 satisfy (jl4j) as 
before (thus forming a maximum IS-multiflow) and take zero values on the auxiliary 
arcs in (Xi), as required. Do similarly for and X3. 

The task of improving the flows within the subgraph (ly) = (ly, A^) is a bit more 
involved. First of all we modify 512 (and 521) so as to get 

(17) gvi{v,v') \ g2z{v,v') ^ gj,x{v,v') = ^ for each f e IV 

(this situation is technically simpler). This is done by decomposing the above- 
mentioned symmetric circulation Aq in (ly) into the sum of an integer circula- 
tion ijj and its symmetric circulation uj' and then by updating 512 := 512 + and 
521 := 521 + (with ijj^uj' extended by zeros to E\ — Aq). Then the equality 
A(f,u') Y,ij 9ijiv,v') = provides dHj). 

The process of improving the flows within (W) consists of 0{W) iterations (the 
idea is borrowed from the algorithm for digraphs in ^U])- At a current iteration, we 
choose a node v E W where some gij{v,v') is nonzero. W.l.o.g, one may assume that 
gi2{v,v') > and gi3{v,v'), g23{v,v') < 0. Let tq := \gi3{v,v')\ and n := \g23{v , v% 
then 5i2(f , v') = tq + vi, by (fTTj) . Let B be the set of (four) auxihary arcs connecting 
V and v'. 

First of all we represent 512 as the sum of two integer Si-Sj flows go,gi such that 
go{v,v') = To and gi{v,v') = ri. To do so, replace B by new terminals t,to,ti and 
arcs a = {t,v'), = (f,to) and ai = (f,ti), and add an arc b from to Si. Define 
h{ai) := Tj, i = 0,1, h{a) := vq + ri, h{h) := val(5i2) and h{e) := 512(e) for the 
remaining arcs e. This turns h into a flow from t to {to,ti}, and we decompose it 
into the sum of integer flows hQ,hi, from t to to and from t to ti, respectively. These 
ho, hi determine the desired 50, 5i in a natural way. 

Combine / := 5o + 5i3- Then f{v,v') = goiv,v') + gi3{v,v') = 0. Update /(e) := 
for each e E B and decompose the updated Si-{s2, S3} flow / into the sum of integer 
flows /o, /i3, from si to S2 and from si to S3, respectively. Then val(/o) = val(5o) and 
val(/i3) = val(5i3). 

Doing similarly for the flow 523 and the flow g[ symmetric to 51 (which have the 
source S2 in common), we obtain corresponding S2-S3 flow /23 and S2-s[ flow f[. 
Finally, update 512 := /0 + /1 (where /i is symmetric to f[), 513 := /13 and 523 := /23- 
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The updated flows gij together with their symmetric ones satisfy gij{e) = for each 
arc e e -B. 

Then we choose a next pair of node mates in W, and so on. Upon termination 
of the process, the resulting flows Qij take zero values on all auxiliary arcs, and it 
is easily seen from the construction that yal{gij) preserves for all pairs ij. So their 
restrictions to E form a maximum IS-multiflow in {G,T,c), as required. 

The above algorithm runs in 0{(j){V, E)) time plus the time needed to perform 
0{W), or 0{V), flow decompositions during the iterative process at Stage 3 (the other 
operations including those in 0(1) symmetric decompositions take 0{VE) time). 
Each of these decompositions is applied to a flow with 0(1) sources and sinks, and 
we use the procedure in Section^lto implement it in 0(£'log(2 + V'^/E)) time. This 
gives the bound 0(V^-Elog(2 + V'^/E)) for the six (or four) terminal case. 

4.2 General case. 

We now describe the algorithm for an arbitrary jS*! > 4. It is based on a recursive 
network partition approach. 

For a current inner Eulerian skew-symmetric network = (G,T, c), with T = 
SUS', the network partition procedure partitions S into two sets 5*1, 5*2 such that l^il = 
[|S|/2] and |^2| = LI'5|/2J and finds a symmetric subset X cV with XnT = SiUS( 
whose induced cut S{X) has minimum capacity c{S{X)). This is done by finding a 
minimum capacity cut S{Y) with y fl T = 5*1 U 5"^ in the underlying undirected 
network for (G, c), and by making the symmetrization X := Y U Y' (relying on 
c{5{Y U Y')) + c{6{Y n Y')) < c{5{Y)) + c{5{Y')) = 2c{6{Y))). 

Next we shrink the subgraph {V — X) of G into two new (extra) terminals ti,t[, 
making each arc in 6°^^{X) enter t[, and each arc in 5™(X) leave ti. Similarly, (X) 
is shrunk into extra terminals ^2,^2? each arc in 5™(X) becomes entering ^2 and each 
arc in 6°'^^{X) becomes leaving ^2- This produces two smaller inner Eulerian networks 
N, = (G, = {V,, Ei),Ti, a) with Ti = SiU S[ U {t„ t^}, ^ = 1, 2, satisfying 

(18) |T,|<4|T|/5, \Vi\<\Vl |^^|<|^|, and + |\/2| = |^| + 4 

(since |Ti| = 8 when |T| = 10). Also for Xi := {tj,t^}, the cut 5{Xi) of Gj has 
minimum capacity among the cuts separating {ti,t[} and Tj — {ti,t[}. 

One application of the network partition procedure, to a current A^, takes one 
minimum cut computation, so it runs in 0{(j){y,E)) time. 

Let Fi be a (recursively found) maximum free IS-multiflow in Ni. The aggregation 
procedure transforms Fi,F2 into a maximum free IS-multiflow F in A^. The flows 
in Fi going from Si to t'i are combined into one (multisource) flow fi from Si to 
t'i, and symmetrically, the flows from ti to S'i are combined into one flow fl. By 
the maximality of Fi and the minimality of Ci{6{Xi)), fi saturates 5™(Xj) and fl 
saturates 6°^*'{Xi). We glue together (the images of) /i and /j, obtaining 5'i-5'2 flow 
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/ in N, and do symmetrically for f2, f[, obtaining /'. These /, /' are decomposed 
symmetrically into a symmetric collection of integer one-source-one-sink flows. Then 
the flows formed from /, /' together with the remaining flows in Fi,F2 connecting 
pairs of terminals in Ti or in T2 give the desired F. (The maximality of F follows 
from the fact that for each s & Si, the total value of flows in Fi leaving s or entering 
s' is equal to the minimum capacity of a cut in A^i separating {s, s'} and Tj — {s, s'}, 
and similarly for 5*2 . The above construction maintains such an equality for F and 
each s E S.) 

At the bottom level (IS*! = 3), we apply the algorithm described in 14.11 

One application of the aggregation procedure, to current A^i, N2, takes 0{SiEi + 

S2E2) time to create the flows /i, /2 as above plus OiVE) time to decompose /, or 

0{VE) time in total (in view of (ITHD). 

It remains to explain that the resulting multiflow F in the initial network N can 
be efficiently transformed into a maximum integer free multiflow in the corresponding 
bidirected network {H,S,'c). We show that 0{VE\ogT) time is sufficient to create 
from F a corresponding sjTiimetric collection {V, a) of weighted T-paths in N; these 
paths determine weighted T-walks in H forming an optimal solution to problem (P) 
with {H,S,c), by the relationship explained in Section IHl We assume that each flow 
/ in F is explicitly given only within its support supp(/) := {e E E : /(e) 7^ 0}. 

Let T be the binary rooted tree formed by all networks arising during the re- 
cursion, with the natural ordering on them. The height of T (or the depth of the 
recursion) is O(logT), in view of the flrst inequality in (fTHj) . For a network N in 
T, let A{N) be the set of terminals from the initial T that are contained in N, and 
F{N) the set of flows in F with both terminals in A{N). We use the fact that for 
incomparable A^', N" in T, the supports of flows in F{N') are disjoint from those in 
F{N"). (Indeed, for the closest common predecessor N of N',N", the minimum cut 
found by the network partition procedure for N separates A{N') and A{N") and is 
saturated by the flows not in F{N') U F{N").) 

We proceed as follows. For each non-leaf network N with children A^i, A^2; combine 
the flows in F with the source in A{Ni) and the sink in A{N2) into one multiterminal 
flow / (in the initial network), and then decompose / into a set V{N) of weighted 
paths from A{Ni) to A{N2) (the circuits appeared in the decomposition are removed). 
This takes 0(V^supp(/)) time. Taken together, the sets V{N), their symmetric sets 
and corresponding paths appeared by decomposing the flows in F having both ter- 
minals in one leaf network, constitute the desired symmetric collection {V,a). To 
estimate the complexity, consider the networks N at height i in T. They are incom- 
parable, so the supports of flows / as above in them are pairwise disjoint. Hence to 
form the sets V{N) for these N takes 0{VE) time in total. This gives the bound 
0{VElogT) for the whole procedure, as declared. 
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4.3 Complexity of the algorithm 

We show that the above algorithm runs in 0{VE i{V, E) logT) time, where £(V, E) : = 
max{l, ln(V"^/-E')}, assuming (f){n,m) = 0(nmlog(2 + n^/m)) (as in Goldberg- 
Tarjan's max flow algorithm). We use induction on the height h{T) of the binary 
tree T (it depends only on |T|). When h{T) = (i.e. |T| = 6), the required time 
bound was shown in 14.11 

Let h{T) > 1 and let Ni,N2 be the children of in T. For i = 1,2, we have 
h{Ti) < h{T) — 1, and by induction the time Xj of the algorithm to solve the problem 
for Ni is bounded from above as 

Ti < Cnimd{ni,mi)\ogTi 

for some appropriately chosen constant C > (specified later). Here rti := \Vi\ and 
TTii := Ei, keeping notation from 14.21 The network partition and aggregation proce- 
dures applied to take time 0{(j){n,m)) and 0{nm), respectively, or Dnmi{n,'m) 
time together, where D is some constant > 0, n := \V\ and m := \E\. Therefore, the 
time r to solve the problem for is estimated as 

(19) T < C{nimii{ni,mi) logTi + n2m2i{n2, ^2) logT2) + Dnmi{n, m). 

We have l{ni,mi) < i{n,mi) (since rii < n) and ■mii{n,mi) < ml{n,m) (this 
follows from rrii < m and from Ina > ln(a6) for 6 > 1 and In a > 1). Also 
ni + n2 = n + 4 and logT^ < logT - log f , by JIHl). Then ^ implies 

/ 5 5\ 

(20) T < C i nm log T — nm log - + 4m log T — 4m log - j i(n,m) + Dnmi{n, m) . 

Since n > \T\ and |r;,m log | grows faster than 4m log n, one can choose constants 
no and C (depending on D) such that the right hand side value in (j20|) becomes 
smaller than Cnmi{n, m) logT for any n > no- (For the networks with \V\ < no, the 
problem is solved in 0{E) time.) This yields the desired time bound. 

Theorem 4.1 A maximum IS-multiflow (resp. a maximum integer free multiflow) 
in an inner Eulerian skew-symmetric (resp. bidirected) network [G = {V,E),T,c) 
can be found in 0{VE\ogT\og{2 + V'^/E)) time. 

5 Fast Flow Decomposition 

For a fixed k G Z+, we consider the problem: 

(D) Given a flow (integer flow) f from S to T, with \S\ + |T| = k, in a digraph 
G = (V, E), find a decomposition f = ^(/^t : s & S,t & T) , where each fst is a 
flow (resp. integer flow) from s to t, 
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and show the following (allowing parallel arcs in G and assuming = 0(E)). 
Theorem 5.1 (D) can be solved in 0{E log{2 + V'^/E)) time. 

Note that when G is acyclic, a decomposition (into one-source-one-sink flows or 
into weighted paths) of any flow in G is carried out in 0{E) time by using a topological 
sorting of the nodes. Sleator and Tarjan JH] showed that any flow / in an arbitrary 
digraph can be decomposed, in 0{E\ogV) time, into a circulation and a flow whose 
support induces an acyclic subgraph of G (so a decomposition of / into one-source- 
one-sink flows can be found with the same complexity 0{E log V)). The algorithm 
in uses sophisticated computational tools, so-called dynamic trees. 

Our approach to solve (D) is based on a node splitting technique and uses only 
simple data structures. Let 11 be the set of pairs st with s E S and t eT. 

In the beginning of the algorithm, we delete from G the arcs e with /(e) = 0. 
Also we sort the nodes v by increasing their degrees deg(f). (This takes 0{E) time.) 
The algorithm applies \V\ iterations. 

At each iteration, we choose a node v with deg(f ) minimum in the current graph 
G = {V, E) . First of all we scan the arcs incident with v to select parallel arcs among 
them. Each tuple of parallel arcs is merged into one arc (and the flows on these 
are added up). The node degrees and the ordering on V are updated accordingly. 
(This preliminary stage is performed in 0{deg{v)) time. As a result, the degree of v 
becomes less than 2|y|.) Then we make at most deg(f) splittings at v. 

More precisely, at a current step of the iteration, we choose an arc entering v and 
an arc leaving v, say, e = {u, v) and e' = {v, w). If e or e' is a loop, we simply delete it 
from G. Otherwise define e := min{/(e), /(e')}. The splitting- off' operation applied to 
(e, V, e') creates a new arc e" from u to w, assigns f{e") := e, updates /(e) := /(e) — e 
and /(e') := f{e') — e, and deletes from G the arc (or arcs) for which the new value 
becomes zero; it takes 0(1) time. The ordering on V is updated accordingly (in 0(1) 
time). Clearly the operation maintains both the divergency at each node and the flow 
integrality (when the original flow is integer). Also div(f) decreases and the number 
of all arcs does not increase. 

At the next step of the iteration, the operation is applied to another pair of arcs, 
one entering and the other leaving v, and so on until such pairs no longer exist. After 
that, if V ^ S UT, then v is removed from G (as div/(f ) = implies deg(w) = 0). 

At the next iteration, we again choose a vertex where the current degree is mini- 
mum, and so on. One can see that after \ V\ iterations, each arc of the resulting graph 
G goes from a source s G 5* to a sink t E T. The decomposition D = {/^t : st G 11} 
for the resulting / in G is trivial: fst{s,t) := f{s,t) and /st(e) := for e 7^ (s,t) 
(letting fst := if the arc {s,t) does not exist in G). 

Now going in the reverse order and applying the corresponding restoration proce- 
dure reverse to the splitting-off one, we transform D into the desired decomposition 
of the initial flow. More precisely, consider a current graph G and the arcs e = {u, v), 
e' = {v,w), e" = {u,w) as above, and let fst, st G 11, be the flows already obtained 
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for the graph G' formed from G by the sphtting-off operation w.r.t. (e, f,e'). For 
each st e 11, add fsti.^") to /st(e) and to fsti.^') and then delete e". (The backward 
iteration concerning v finishes with restoring the corresponding tuples of parallel arcs 
incident with v and assigning, in a due way, the flows fst on these arcs.) Eventually, 
we obtain the desired decomposition {fst : st G 11} of the initial /. (Strictly speaking, 
we have g := J2stfst < / a^id div f_g{v) = for all v E V; so one should add the 
circulation f — g to one of the flows fst-) 

Next we estimate complexity of the above algorithm. Let Vi,V2,... ,v^v\ be the 
sequence of nodes in the splitting-off process. Since |n| = 0(1), the restoration 
process is only 0(1) times slower than the splitting-off one. (This is just where 
we essentially use the condition that / has 0(1) terminals.) Using this fact, one 
can conclude that the algorithm runs in 0{E + A) time for the initial E, where 
A := deg*(f i) + . . . + deg*{v\v\) and deg*(T;) denotes the degree of v at the beginning 
of splitting at v. Each iteration i in the former process does not increase the number of 
arcs of the current graph and decreases the number of nodes by one, unless Vi G SUT. 
So deg*(fj+i) is at most 2|i?|/(|l^| — i). Summing up the latter numbers over i, we 
obtain A = 0{E log V), which is worse than the time bound in Theorem 15.11 

However, we can estimate A more carefully, by using the inequality deg*(fj) < 
2{\V\—i + k + l) (provided by merging parallel arcs incident with Vi). For any integer 
1 < A < \V\, apply the first bound on deg*(fj) for i = 1, . . . , |V^| — A, and the second 
bound for 2 = |V^| — A + 1, . . . , \V\. This gives 

or A = 0{E \og{V/X) + A^). Now taking A := min{|V|, [v^D' obtain A = 
0(Elog(2 + V'^/E)), and the theorem follows. 

6 Fast Skew-Symmetric Flow Decomposition 

In this section Theorem l5.1l is extended to (skew-)symmetric flows. For a fixed k G Z+, 
we consider the problem: 

(DS) Given an integer symmetric flow f from S = {si, . . . , Sk} to S' = cr{S) in a 
skew-symmetric graph G = {V, E), find a decomposition of f of the form 

where each fij is an integer flow from Si to s'j and f^j is symmetric to fij . 

Note that f-j is a flow from Sj to s'^. So in the above decomposition, for i,j G 
{1, . . . , k}, Si and s'j are connected by the only flow fij ii i < j, by only /jj if i > j, 
and by the two flows fa and f■^ if i = j. We show the following 
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Theorem 6.1 (DS) can he solved in 0{E\og{2 + V'^ / E)) time. 

This generahzes Theorem lS.ll for integer flows because a digraph D with an integer 
S-T flow g is turned into a skew-symmetric graph with an integer symmetric (SUT')- 
(5" U T) flow by adding a disjoint copy of the reverse to D with the flow reverse to g 
in it. 

Our algorithm to solve (DS) relies on the following lemma (where, as before, 
primes are used for the corresponding mate objects). 

Lemma 6.2 Let g he a (not necessarily symmetric) half-integer flow from S to T in 
a skew-symmetric graph G = {V,E) such that divg{v) is an integer for each v & V. 
Let g + g' he integer. Then there exists, and can he found in 0{E) time, an integer 
flow h in G such that h + h' = g + g' and div/j(t') = divg(f ) for all v ^ V. 

Proof. Let Eq be the set of arcs e with g{e) ^ Z. The integrality of g -\- g' implies 
E'q = Eq, so the subgraph (Eq) induced by Eq is skew-symmetric. Also the half- 
integrality of g and the integrality of div^ imply that each node is incident with an 
even number of arcs in Eq. So the underlying undirected graph H of (Eq) is Eulerian. 

We grow a (simple) path P in (Eq) such that PCiP' = 0, starting with an arbitrary 
node Vq and allowing backward arcs in P. Let v be the last node of the current P, 
and choose an arc e & Eq incident with v and different from the last arc of P (e exists 
as H is Eulerian). Let u be the end of e different from v. Three cases are possible, (i) 
If both u, u' are not in P, we increase P by adding e, u, and continue the process, (ii) 
If G P, we remove the part of P from u to f , obtaining the new current path from 
fo to M, and add e to the removed part, forming circuit C (with possible backward 
arcs), (iii) If u' G P, we remove the part Q of P from u' to v, obtaining the new 
current path, and add e, Q' and e' to Q, forming circuit C (which is reverse to G'). 

In case (ii), we update g by pushing half-unit along G (i.e., by setting g{a) := 
5f(a) + I for the forward arcs a in G, and g{a) := g{a) — ^ for the backward arcs a) 
and by pushing half-unit along the circuit reverse to G'. Accordingly, we update g' 
by pushing half-unit along G' and along the circuit reverse to G. And in case (iii), 
g {g') is updated by pushing half-unit along G (resp. G'). In both cases, the new g' 
is symmetric to the new g and each of the functions g + g' and div^ preserves. Also 
Eq decreases by the set of arcs occurring in C U G', and the new H is Eulerian. We 
continue the process with the new P. 

The final g,g' give the desired h, h' . The bound 0{E) is obvious. ■ 

Return to problem (DS). Add to G new nodes t,t' and arcs (t, s,) and (s^,t'), 
forming skew-symmetric graph Gi, and extend / to an IS-flow from t to t' in Gi in a 
natural way. The fact that / is integer and symmetric implies that divj(t) is even. 

So we can apply Lemma W7]\ to the flows g := g' := obtaining corresponding 
integer flows /i, h' . The restriction h of h to E is an integer flow from S to S", and we 
apply the 0{Elog{2 -\- l^^/£^))-algorithm from Section Elto decompose it as 
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where hij is an integer flow from Sj to s'j. Then the flows fa := ha for i = 1, . . . ,k, 
and fij := hij + h'-^ ioT 1 < i < j < k required, and Tlieorem Ifi. II follows. 

Remark. The above proof involves the following corollary from Lemma 16.21 

Corollary 6.3 Let f be an IS-flow from S to S' in a skew-symmetric graph G = 
{V,E) (where \S\ is not fixed), and let div/(f) be even for all v E V. Then there 
exists, and can be found in 0{E) time, an integer flow g from S to S' such that 
f = 9 + 9' (^''^d div/(t>) = 2divg(f ) for all v E V . 
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